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FOREWORD 


Since  the  early  1960s  artificial  earth  satellites  have  been  tracked 
using  the  Doppler  principle.  Until  the  early  1970s  a  large  part  of  this 
data  was  processed  in  orbit  determination  programs  as  frequency  data,  a  pro¬ 
cedure  which  resulted  in  a  small  error  due  to  a  truncated  time  series.  Today, 
the  data  fran  only  a  few  stations  is  processed  by  this  method.  This  report 
describes  a  better  way  to  process  this  data  and  gives  the  differences  in 
results  for  the  two  methods. 

This  work  was  accomplished  in  the  Satellite  Geodesy  Branch,  Astronautics 
and  Geodesy  Division.  The  report  was  reviewed  by  R.  J.  Anderle,  Head,  Astronautics 
and  Geodesy  Division. 
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INTRODUCTION 


Doppler  observations  of  artificial  earth  satellites  have  been  taken  since 
the  early  1960s  to  determine  the  satellite's  orbit  if  the  positions  of  the 
observing  sites  are  known  or  to  determine  the  position  of  the  observing  site 
if  the  satellite  orbit  is  known.  The  Doppler  observations  are  made  at  the  ob¬ 
serving  stations  by  the  station  receiving  a  continuous  wave  signal  generated 
by  a  satellite  at  a  known  frequency.  1116  received  frequency  is  beat  against 
a  frequency  generated  by  a  precise  ground-station  oscillator.  The  method  of 
counting  these  beat  cycles  differs.  Until  1971  most  of  the  data  was  taken 
by  counting  a  preset  number  of  cycles  and  recording  the  time  required  for  the 
count.  The  count  was  typically  done  once  every  4  seconds,  and  the  number  of 
cycles  to  count  was  chosen  so  that  the  count  interval  would  last  just  under 
1  second.  This  type  of  data  is  usually  labeled  as  Sampled  Doppler.  Today 
only  the  OPNET  stations  located  in  Maine,  Minnesota,  California,  and  Hawaii 
take  data  using  this  method,  for  in  1970-71  the  other  stations  were  converted 
to  a  more  automatic  mode  of  operation.  As  part  of  this  change,  the  counting 
device  was  modified  so  that  the  time  for  the  completion  of  one  count  is  now 
identical  to  the  time  for  the  start  of  the  next  count.  The  number  of  cycles 
to  be  counted  is  set  so  that  the  counting  interval  lasts  from  10  to  30  seconds. 
This  type  of  data  is  labeled  Continuous-Count  Integral  Doppler  (CCID).  The 
Sampled  Doppler  data  has  been  processed  as  frequency  data  in  orbit  determination 
programs,  although  the  technique  for  doing  this  leads  to  an  error  in  position 
due  to  an  approximation  when  the  number  of  cycles  counted  in  an  interval 
divided  by  the  time  interval  is  treated  as  a  frequency  in  the  equations  of 
condition.  This  truncated  series  results  in  an  error  that  has  previously 
been  neglected.  The  behavior  of  this  error  was  evaluated  by  treating  this 
Sampled  Doppler  data  as  range  difference  rather  than  frequency.  This  report 
gives  the  magnitude  of  the  error  involved  and  a  way  of  avoiding  it. 


PROCEDURE 


Filtered  frequency  observations  and  a  readily  modifiable  station  position 
determination  program  were  initially  used  for  this  study.  Their  use  provided 
the  most  expedient  way  to  obtain  the  magnitude  of  the  error  to  determine  if 
it  was  of  enough  consequence  to  pursue  and  to  justify  modifying  our  large- 
scale  orbit  determination  program,  CELEST.*  Here  we  had  available  the  time 
of  the  midpoint  of  the  interval  over  which  the  Sampled  Doppler  data  was  taken, 
and  a  corresponding  frequency  observation  that  had  been  computed  from  the 
ratio  of  the  Doppler  count  to  the  time  interval.  Needed  were  the  times  at 
the  beginning  and  at  the  end  of  the  counting  interval  and  a  way  of  converting 
the  frequency  value  of  observation  to  a  range  difference  observation. 


*  James  W.  O'Toole,  CELEST  Computer  Program  for  Computing  Satellite  Orbits , 
Naval  Surface  Weapons  Center  Technical  Report  TR-3565,  Dahlgren,  Virginia, 
October  1976. 
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To  obtain  the  length  of  the  count  interval 


<fR  -  fS>  -  Af 


where  At  *  time  to  count  Nc  cycles 

Nc  =  number  of  cycles  to  count 
fR  =  station  reference  frequency 
fg  =  transmitted  frequency 

Nc 

Af  =  frequency  data  observation  value  =  (fR  -  — ) 


f 


S 


from  which  the  begin  and  end  times  of  the  sample  interval  are  obtained. 
The  range  difference  value  may  then  be  computed  by 


AR  =  7  [Nc  -  <fR  -  fS>  At] 

§ 

where  AR  =  range  difference  observation  value 
c  =  speed  of  light 

This  method  was  sufficient  to  indicate  that  there  was  an  error  of  signifi¬ 
cant  magnitude  to  justify  the  modification  of  our  main  orbit  computation  program, 
CELEST.  The  best  way  to  do  the  conversion  is  by  using  the  raw  observational 
data  because  one  has  available  the  actual  time  required  to  count  the  preset 
number  of  cycles  and  the  time  the  count  started.  This  is  the  way  the  conversion 
will  be  made  in  practice  and  the  method  used  for  the  later  parts  of  the  study. 

The  interval  over  which  the  sample  was  taken  was  varied  to  see  if  the  re¬ 
siduals  of  fit  would  vary.  It  was  expected  that  the  residual  differences  would 
vary  quadratically  with  the  interval  change  since  powers  of  the  reciprocal 
square  of  the  time  interval  and  higher  are  ignored  in  the  frequency  computation. 
The  interval  was  varied  as  the  frequency  data  was  converted  to  range  difference. 
The  range  difference  data  was  then  refit  to  obtain  new  residuals.  The  sample 
intervals  used  were  the  full  interval  of  the  original  data,  one-half  the  in¬ 
terval,  and  one-fifth  the  full  interval. 

Two  days  of  data  from  the  Navy  navigation  satellite,  1970-67A,  were  used 
for  the  test.  This  is  a  polar  satellite  with  a  perigee  height  of  about  970  km 
and  an  eccentricity  of  0.017.  Other  spans  of  data  were  also  used  as  checks 
to  be  sure  the  discrepancies  were  not  span-dependent. 
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RESULTS 


When  changing  from  frequency  to  range-difference  data  the  resultant 
station  positions  varied  by  about  1.4  meters  in  height  with  small  changes 
in  latitude  and  longitude  for  a  sample  interval  close  to  1  second.  When 
the  sample  interval  was  decreased,  the  station  height  residual  differences 
decreased  approximately  as  the  square  of  the  interval  change.  This  change  is 
seen  in  the  pass-by-pass  differences  in  Table  1,  or  more  readily  in  Table  2, 
the  combined  pass  solution  with  the  differences  given  in  three  components. 

These  changes  are  also  seen  in  Figures  1-6  where  the  first  three  show  the 
longitude  differences  as  the  interval  decreases  and  the  second  three  show  the 
corresponding  latitude  changes.  The  longitude  difference  plots  also  show 
on  which  side  of  the  station  the  satellite  was,  and  the  latitude  difference 
plots  show  the  direction  of  the  pass.  These  are  given  to  show  that  there  is 
not  a  direct  correlation  between  the  pass  geometry  and  the  station  position 
differences.  Table  3  and  Figures  7  and  8  give  the  station  navigation  dif¬ 
ferences  in  the  alongtrack  and  radial  frames.  Figure  7  radial  differences 
show  the  bias  caused  by  the  difference  between  processing  as  frequency  and 
range  difference.  Table  3  and  the  corresponding  figures  (7  and  8)  were  produced 
by  CELEST.  CELEST  was  modified  after  the  station  position  differences  were 
shown  to  be  of  a  significant  magnitude. 


CONCLUSION 


The  major  resultant  differences  between  the  two  processing  techniques 
is  that  processing  the  data  as  frequency  produces  station  heights  that  are 
about  1.4  meters  too  low  for  Doppler  count  intervals  of  about  1  second.  In 
producing  today's  satellite  orbits  or  point  positions  that  are  accurate  to 
a  few  meters,  this  error  is  too  large  to  ignore  even  though  there  are  currently 
only  a  few  observing  stations  that  record  this  type  of  data.  In  processing 
data  for  periods  prior  to  1971,  when  most  stations  were  taking  frequency 
data  (such  as  when  doing  a  general  geodetic  solution) ,  this  discrepancy  should 
not  be  overlooked  as  an  error  source. 


LIMITATIONS  OF  THE  METHOD 


Least-squares  solutions  for  geodetic  parameters  involve  the  partial 
derivatives  of  the  data  with  respect  to  the  parameters.  For  range-difference 
data,  partial  derivatives  for  range  are  usually  formed  at  the  beginning  and 
end  of  the  count  interval  and  are  differenced  to  obtain  the  partial  deriva¬ 
tive  for  range  difference.  For  count  intervals  of  1  second  or  less,  computer 
word  length  begins  to  introduce  significant  errors  in  the  partial  derivatives. 
An  approach  taken  to  avoid  this  problem  will  be  the  subject  of  a  later  report. 


3 


Table  1.  Station  Position  Residual  Differences 


Frequency  -  Range  Difference  (a) 


Station* 

Pay 

Hour 

Elevation 

(Deg) 

sir 

Sasple  Interval 

Interval/2 

Interval/5 

321 

60 

1 

63 

Long 

3.55 

1.16 

0.55 

Lat 

-0.42 

-0.16 

-0.10 

330 

1 

25 

Long 

-0.70 

-0.45 

-0.23 

Lat 

-0.17 

-0.06 

-0.01 

311 

11 

64 

Long 

3.58 

1.08 

0.55 

Lat 

-0.06 

-0.01 

0.00 

321 

11 

28 

Long 

-1.52 

-0.53 

-0.30 

Lat 

-0.37 

-0.12 

-0.11 

330 

13 

23 

Long 

-1.38 

-0.49 

-0.27 

Lat 

-0.21 

-0.08 

-0.11 

320 

13 

57 

Long 

2.73 

0.95 

0.39 

Lat 

0.10 

-0.03 

-0.04 

330 

15 

52 

Long 

2.43 

0.82 

0.38 

Lat 

-0.04 

0.00 

0.00 

311 

22 

63 

Long 

-3.22 

-1.01 

-0.44 

Lat 

0.00 

-0.01 

0.01 

320 

61 

0 

70 

Long 

-4.05 

-1.41 

-0.S1 

Lat 

0.16 

0.02 

0.01 

311 

0 

27 

Long 

1.45 

0.50 

0.21 

Lat 

-0.61 

-0.28 

-0.11 

311 

10 

66 

Long 

-4.43 

-1.39 

-0.53 

Lat 

-0.61 

-0.16 

-0.03 

320 

12 

74 

Long 

-6.80 

-2.18 

-0.85 

Lat 

-0.79 

-0.25 

-0.08 

311 

12 

24 

Long 

1.21 

0.41 

0.34 

Lat 

-0.07 

-0.05 

0.02 

311 

22 

32 

Long 

-1.38 

-0.50 

-0.25 

Lat 

-0.20 

-0.08 

0.02 

311 

23 

53 

Long 

2.22 

0.49 

0.04 

Lat 

1.65 

1.89 

1.96 

320 

23 

36 

Long 

-0.44 

0.53 

0.72 

Lat 

1.48 

1.58 

1.61 

*  Station  Locations:  311  Prospect  Harbor,  Maine 

320  Roseaont,  Minnesota 

321  Roaenont,  Minnesota 
330  Laguna  Peak,  California 
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Table  2.  Combined-Pass  Station  Position  Differences 


Frequency  -  Range  Difference  (m) 


No. 

Full 

Station* 

of  Passes 

Sample  Interval 

Interval  /2 

Interval  /5 

311 

7 

Long 

0.30 

0.26 

0.29 

Lat 

0.71 

0.68 

0.68 

Ht 

-1.37 

-0.38 

-0.12 

320 

4 

Long 

0.43 

0.42 

0.36 

Lat 

0.30 

0.28 

0.27 

Ht 

-1.32 

-0.39 

-0.12 

321 

2 

Long 

-0.25 

-0.11 

-0.09 

Lat 

0.12 

0.03 

-0.01 

Ht 

-1.48 

-0.50 

-0.26 

330 

3 

Long 

-0.06 

o 

• 

o 

1 

-0.05 

Lat 

0.00 

0.01 

0.01 

Ht 

-1.52 

-0.53 

-0.25 

*  See  Table  1 
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Figure  2.  Longitudinal  Differences  in  Station  Position  Residuals 
Sample  Interval/2  (Frequency  -  Range  Differences) 
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Figure  3.  Longitudinal  Differences  in  Station  Position  Residuals 
Sample  Interval/5  (Frequency  -  Range  Differences) 
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Figure  4.  Latitudinal  Differences  in  Station  Position  Residuals 
Full  Sample  Interval  (Frequency  -  Range  Differences) 
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Figure  5.  Latitudinal  Differences  in  Station  Position  Residuals 
Sample  Interval/2  (Frequency  -  Range  Differences) 
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Figure  6.  Latitudinal  Differences  in  Station  Position  Residuals 
Sample  Interval/5  (Frequency  -  Range  Differences) 


Table  3.  CELEST  Station  Navigation  Differences 
(Frequency  -  Range  Differences)  (m) 

Frequency  -  Range 
Elevation  Difference  (m) 


Station* 

Day 

Hour 

(Deg) 

Radial 

Along-Track 

321 

60 

1 

63 

-1.2 

-0.2 

330 

1 

25 

-1.4 

0.3 

311 

11 

64 

-0.9 

0.1 

321 

11 

28 

-0.4 

-3.3 

330 

13 

23 

-0.1 

1.9 

320 

13 

57 

-0.3 

-0.3 

330 

15 

52 

-1.6 

0.0 

311 

22 

63 

-1.7 

2.5 

320 

61 

0 

70 

-2.1 

0.3 

311 

0 

27 

4.1 

14.1 

311 

10 

66 

-1.6 

0.2 

320 

12 

74 

-1.6 

-0.4 

311 

12 

24 

-2.7 

1.6 

311 

22 

32 

0.5 

-1.6 

311 

23 

53 

-2.1 

-0.2 

320 

23 

36 

-0.7 

0.1 

*  See  Table  1. 
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